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Abstract: We discuss the comparisons of several kinematic distributions, measured by 
the ATLAS and CMS Cohaborations for dijet production, with the results computed with 
three approaches: using a fixed next-to-leading order calculation, with the POWHEG BOX and 
with HEJ. We study other interesting kinematic distributions, exploring the effect of several 
sets of cuts. Previous experimental studies have dealt with kinematic distributions where 
the predictions of the three approaches were very similar. In this article, we investigate 
kinematic distributions where the resummed effects in POWHEG and HEJ are enhanced with 
respect to the fixed next-to- leading order result, since different QCD-radiation regimes are 
probed, and we suggest a few kinematic distributions where the predictions from POWHEG 
and HEJ are different, and could be further investigated experimentally. 
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1. Introduction 

Dijet production is one of the cornerstone processes at the LHC. The cross section for jet 
production is very large, making it an important testing ground for our understanding of 
QCD at high-energy scales. In addition, jet production is an important background for 
many searches for new physics. It is therefore essential to probe and test our theoretical 
predictions. 

Dijet production is not only important by itself but also when the jets are accompa- 
nying other particles, such as the W jZ vector bosons or the Higgs boson. For example, 
Higgs boson production plus two jets is an important process in the Standard Model (SM). 
Two different processes have Hjj as final-state signature: in the first one, the Higgs boson 
couples to gluons via a top-quark loop [1, 2], while in the second one, the Higgs boson 
is produced by vector-boson fusion (VBF) [3, 4, 5], stemming from two incoming quarks. 
The vector-boson fusion process is particularly important to extract the HWW and HZZ 
couplings, and to confirm if the detected Higgs boson is the scalar responsible for the 
electroweak symmetry breaking. In addition, information on the CP properties of the 
Higgs boson can be extracted by studying the azimuthal distributions of the two hardest 
jets [6, 7], in events with a large (> 3 units) rapidity separation between them. In fact, 
the signature of the VBF Higgs boson production are two jets well separated in rapidity. 
In addition, we expect very low hadronic activity between the two hardest jets, due to the 
exchange of the colorless vector bosons in the t channel [8], on the contrary of what is 
expected for the gluon-fusion production mechanism. A key feature is then the study of 
the efficiency of the central-jet veto, to suppress gluon-fusion processes [1, 2], where it is 
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well known that the higher-order corrections [9, 10, 11, 12] are very significant [13], and a 
detailed understanding [14, 15, 16, 17] of the structure of the radiation pattern is needed. 

Some of the features in Hjj production in gluon fusion are in fact universal [13] to 
dijet processes in general, like W/Z + 2 jets or pure dijet production, and may help as a 
testing ground for techniques which can be applied in the Higgs boson searches and studies. 
Before the possible study of processes with a Higgs boson, it thus becomes interesting to 
investigate the (hard) radiation pattern in events with at least two jets, in particular for 
events with a non-negligible rapidity separation between the two jets. The central question 
is whether a framework based on a (possibly next-to-lcading-order-matched) parton shower 
(which resums the radiation resulting from a large ratio in transverse scales) is sufficient 
for the description of additional jets, or whether BFKL-type [18, 19, 20, 21] effects from 
hard, wide-angle emissions, have already become important at the center of mass energy 
of the LHC (currently 7 TeV). 

In light of this, there have been a number of very interesting experimental studies in 
dijet production by both the ATLAS [22, 23, 24] and CMS [25, 26, 27] Collaborations. It is 
already clear that higher-order QCD contributions, beyond a fixed order, low multiplicity 
calculation, can be important, because the large available phase space for jet emission at 
the LHC compensates for the suppression of extra powers in the strong coupling constant. 

In this paper, we compare two theoretical approaches to dijet production, that include 
higher order effects: POWHEG [28, 29, 30] and HEJ [31, 32, 33]. The POWHEG method suc- 
cessfully merges a fixed next-to-leading order (NLO) calculation with a parton shower pro- 
gram, that resums leading logarithmic contributions from collincar emissions. In this study, 
the POWHEG results obtained with the POWHEG BOX [34] are interfaced with the transverse- 
niomcntum-ordcrcd shower provided by PYTHIA 6.4.25 [35]. In contrast, the starting point 
for HEJ (High Energy Jets) is an all-order approximation to the hard scattering matrix 
element in the regime of wide-angle QCD emissions. HEJ is accurate at leading logarithmic 
precision in the invariant mass of any two jets. This is then supplemented with the miss- 
ing contributions (through a merging and reweighting procedure) necessary to also ensure 
tree-level accuracy for final states with up to four jets. The tree-level matrix elements are 
taken from Standalone Madgraph [36]. 

The layout of the paper is the following: in section 2 we discuss in detail existing 
experimental analyses and theoretical predictions. Since our goal is to compare the POWHEG 
and HEJ approaches to the fixed NLO results too, in section 3 we study several sets of cuts, 
that do not induce a large hierarchy in the transverse scales but impose a large rapidity 
separation between the jets, and we investigate for which set the NLO results remain 
reliable and do not suffer from large unresummed logarithmic terms, due to incomplete 
cancellation of virtual contributions [37, 38, 39, 30]. In section 4, we investigate and 
propose a few analyses which will better expose the differences between POWHEG and HEJ, 
and will allow us to better understand which is the dominant mechanism in the generation 
of radiation. Finally, we summarize our findings in section 5. 
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Figure 1: Both plots from the ATLAS study [22] on the gap fraction, defined as the fraction of 
events with no additional jets in the rapidity region between the tagging jets, as a function of the 
difference of the rapidity of the two jets. In the left-hand side plot, the tagging jets are the most 
forward and most backward jet, while in the right-hand side plot, they are the two hardest jets in 
each event. In both plots, pT is the average transverse momentum of two tagging jets. All jets are 
required to have pj. > 20 GeV and absolute rapidity \y\ < 4.4. 



2. Dijet production at the LHC 

In this section, we discuss in detail existing experimental analyses where dijet data collected 
by ATLAS and CMS have been compared with the theoretical results obtained with the 
POWHEG BOX and HEJ. The POWHEG and HEJ approaches are clearly very different in their 
description of QCD radiation: one approach resums collinear emissions while the other soft 
and hard, wide-angle emissions. Nevertheless, for several kinematic distributions (see for 
example ref. [22]) the predictions from POWHEG and HEJ are very similar. This is due, as 
we show in this section, to the inclusiveness of the studied kinematic distributions and on 
the values of the cuts applied in the experimental analysis. 

2.1 ATLAS results 

Among several dijet results published by the ATLAS Collaboration, one is of particular 
interest to us: a study of the production of additional jets from a dijet system [22]. In 
this study, jets are reconstructed using the anti-A;T jet algorithm [40] with R = 0.6 and 
required to have a transverse momentum above 20 GeV, with absolute rapidity less than 
4.4. We show in fig. 1 extracts from this ATLAS study, where the gap fraction, defined as 
the fraction of events with no additional jets in the rapidity region between the two tagging 
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jets, is plotted as a function of Ay, the difference of the rapidities of two tagging jets. In 
the left plot, the tagging jets are defined to be the most forward and most backward jet 
(in rapidity), while, in the right plot, they are chosen to be the two hardest jets (highest 
transverse momentum) in each event. In both plots, pT is the average transverse momentum 
of two tagging jets, and results are shown for slices of the average transverse momentum 
Pj, ranging from 70 GeV up to 500 GeV. The experimental data is then compared with 
the predictions from HEJ (the blue band indicating the scale variation obtained by varying 
the (equal) renormalization and factorization scale by a factor of two), and the POWHEG 
predictions showered by PYTHIA and HERWIG. 

We would like to add a few comments to these findings: 

1. In the left plot, the theoretical predictions obtained in the two very different ap- 
proaches of HEJ and POWHEG are very similar, and agree with data over a wide range 
of rapidity intervals and average transverse-momentum slices. For this kinematic 
quantity, where the tagging jets are selected to be the forward/backward ones, a 
large hierarchy in the transverse momenta of these jets develops as the average trans- 
verse momentum of the forward/backward jets is increased, so that at least one of 
the forward/backward jets must be very hard. Indeed, it is observed that for this se- 
lection, the average difference in the transverse momentum of the forward/backward 
jet increases systematically with increasing Ay. Since the veto scale {Qo in the label 
of the figure) for counting additional jets is much smaller than the average transverse 
momentum of the tagging jets, the jet production is driven by relatively-soft emis- 
sion from the dijet system. Both HEJ and POWHEG, with different approximations, do 
include these multiple emissions giving a good agreement with data. 

2. In the right plot of fig. 1, where the tagging jets are chosen to be the two hardest jets 
in the event, as the average of the two hardest jets increases to 5 times or more of 
the veto scale Qq, the HEJ prediction starts deviating from data, underestimating the 
amount of radiation (i.e. the prediction for the gap fraction is larger than the data). 
This behavior is expected, since the component of events added with naive tree- 
level matching increases with increasing pr- This component receives no systematic 
treatment of soft resummation within HEJ, a situation which would be improved by 
a complete matching with a parton shower. Progress in this direction has recently 
been made in ref. [41]. 

The POWHEG description includes the effects of collinear emissions through the shower 
formulations, and the theoretical predictions perform well for both the two kinematic 
distributions in fig. 1 (particularly when using the PYTHIA shower). However, as can 
be seen in the right plot, for larger rapidity spans and modest pT, the POWHEG de- 
scription undershoots the data. Indeed, POWHEG contains no systematic resummation 
of all the leading logarithmic terms for large Ay. Overall, the study reported by 
ATLAS shows best agreement with the predictions of POWHEG+PYTHIA, but all the 
studies involve a hierarchy of transverse scales, and therefore, by construction, will 
favor the description with the systematic collinear resummation of a parton shower. 
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Note that the results for PDWHEG+HERWIG are consistently below data (i.e. the events 
contain too many jets). The differences between the results from POWHEG+PYTHIA and 
POWHEG+HERWIG should be considered as a theoretical uncertainty connected to the 
different shower algorithm. 

As a final comment, we are interested in separating the two mechanisms for creating ad- 
ditional jet activity, i.e. a hierarchy of transverse scales (as in the POWHEG shower), and 
the opening of phase space as the rapidity span between two jets is increased (as in HE J). 
For this separation to take place, it is obviously necessary to introduce a selection criteria 
which does not automatically generate a hierarchy in the transverse scales as the rapidity 
span increases. This is partially done in the study performed by the CMS Collaboration, 
discussed in the following section. 

2.2 CMS results 

CMS has reported a study [27] on dijet production with just a simple selection criteria on 
the transverse momenta of jets. Jets are reconstructed using the anti-fcx algorithm with R = 
0.5, and are requested to have Pt > 35 GeV. Events are then required to contain at least 
one forward jet (3.2 < |77-^| < 4.7) and at least one central jet {\r]'^\ < 2.8), where r) is the 
pseudorapidity of the jets. The transverse-momentum spectrum of the hardest central and 
hardest forward jet is then studied. Obviously, any difference in the two spectra is a result 
of radiation beyond the trcc-lcvcl description of back-to-back partons. While the CMS 
study extends the pseudo-rapidity region of jets up to 4.7 units, the transverse momentum 
distributions are integrated over these pseudo-rapidity ranges. Crucially, however, no large 
Px-hierarchy is induced by the cuts, and this gives a cleaner study of the separate effects 
of the relatively modest rapidity gap. 

In the CMS analysis, HEJ describes the px-spcctrum well for both the central and 
the forward jets. POWHEG+HERWIG describes the shape correctly, but the normalization 
is consistently high. The POWHEG+PYTHIA description of the forward jet px-distribution 
performs well, but the description for the central jet shows deviations in both shape and 
normalization. As the events in this analysis have been specifically selected to have a 
non-negligible rapidity span, this slight deviation could be attributed to the absence of a 
systematic treatment of the dominant logarithmic terms for increasing A77. 

2.3 Summary 

The analyses discussed so far show that the descriptions of both POWHEG and HEJ are 
performing well, and in broad agreement. The close agreement between the two can, to 
some extent, be attributed to the requirement of a large px-hierarchy in the study by 
ATLAS, or the modest average rapidity spans in the study by CMS. 

In the rest of this paper, we investigate various observables which can expose the 
differences among the fixed NLO calculation, the POWHEG and HEJ approaches. The first 
task is therefore to develop a set of cuts, for which the NLO prediction for dijet production 
is physically meaningful. This is the topic of the next section. 
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Figure 2: The transverse momentum distributions of the leading forward jet (3.2 < \ri^ < 4.7) 
and the leading central jet {\rf\ < 2.8) in a sample which requires at least one jet with > 35 GeV 
in each region. The top and bottom row contain the same data points, but different theoretical 
predictions. The plots are taken from [27]. 



3. Reliability of the NLO predictions 

It has been known for a while that fixed-order results for dijet production are not reliable 
when symmetric cuts are applied to the transverse momentum of the two hardest jets. This 
was first noticed in [42, 37], in the context of electron-proton collisions at HERA, and later 
also in hadron collisions [38] . A detailed theoretical discussion of the origin of this fact can 
be found in [39], where the next-to-leading logarithmic resummation of soft logarithms was 
also performed. 

While the fixed-order theoretical predictions display an unphysical behavior as the 
symmetric-cut limit is reached, the experimental data are obviously not affected by that. 
If we are interested in a comparison of the fixed NLO result with data, then we have to 
impose asymmetric cuts 

P't > Pt" , Pt > Pt'" + Apt , Apt > , (3.1) 
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Figure 3: The dependence of the da/dAyn, distribution on the asymmetry of jet cuts (left plot) 
and of the total cross section in different Aytb slices, as a function of the jet cuts asymmetry Ap-r 
(right plot). 



i.e. all jets are required to have a minimum transverse momentum p™™, while a stronger 
constraint is applied to the hardest jet in the event, with transverse momentum p^^. In 
our notation, Apx quantifies then the asymmetry on the cuts. Since in this paper we are 
interested in studying the size of QCD corrections in regions traditionally used to probe 
BFKL-like effects, we allow for quite forward jets. Therefore we impose the condition, 
lUjl < 4.7, on the jet rapidities, and we set p™™ = 35 GeV (similar to the CMS cut 
discussed in the previous section). Jets are reconstructed using the anti-Zcx jet algorithm [40] 
with R = 0.5 and default recombination scheme. The NLO results, as well as the POWHEG 
distributions shown in section 4, have been obtained with renormalization and factorization 
scales set to the POWHEG underlying-Born transverse momentum, i.e. the px of the partons 
in the Born-like 2 — )• 2 kinematics, starting point for the generation of radiation [30]. For 
all the plots in the paper we have used the MSTW2008 [43] parton distribution function 
set. The aim of this section is to explore several sets of cuts to be applied to the jets and 
to clarify the breaking-point, where the NLO calculation becomes unreliable. 

In fig. 3 we illustrate the effects of the cuts of eq. (3.1): in the left plot, we display the 
differential cross section for dijet inclusive production as a function of Ayf^, the difference 
in rapidity of the most forward and most backward jet. The various curves correspond to 
different values of Ap^ in eq. (3.1), ranging from Ap^ = 1 GeV up to Ap^ = 65 GeV. On 
physical grounds, we expect a decrease of the differential cross section as Ap^ increases. 
However, and as noted in earlier studies [30], we observe that the results do not show this 
behaviour. Indeed, over a wide range in Ayt^,, the cross sections increase for increasing ApT 
till a maximum is reached at ApT ~ 5 GeV. Then, for any further increase of Ap^; the cross 
section decreases, as expected. This is caused by large logarithmic term in ApT, arising 
from a veto in the emission of radiation above the Ap^ scale, which causes uncanceled 
virtual corrections to build up above this scale. The same behaviour is also evident in the 
right plot of fig. 3, where the cross section is plotted against ApT for various slices in Aj/f^. 
The presence of uncanceled virtual corrections becomes manifest for small Apx, where the 
cross section is unphysically negative. 
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Figure 4: The differential cross sections da/dAyn, for several values of p* at X = 2/3 (left plot) 
and for several values of X for p* — 75 GeV (right plot). For both plots, jets have transverse 
momentum greater than 35 GeV. 

In fig. 4, we plot the differential cross section as a function of Ayj^ using another set 
of cuts, BFKL inspired: 

pi,>pf\ P^'+P^>p*^, ^>X, |2/,|<4.7, (3.2) 

Pt 

where ji and j2 denote the hardest and next-to-hardest jet. In the left plot of fig. 4, we 
show the cross sections obtained at a fixed value of X (chosen here to be 2/3) for several 
values of p* ^. The role of Ap^ in fig. 3 is now played by the distance — 2p™™. In 
fact, when p^ reaches its minimum value equal to 2p™™, jets can approach the symmetric- 
cut configuration, exposing again large logarithms in a fixed NLO calculation. The red 
curve shows precisely the unphysical behavior when p^ = 2p™™ = 70 GeV, giving rise 
to a negative cross section. We have checked that this conclusion holds regardless of the 
value of X chosen. As p^ increases beyond its minimal value, the cross sections exhibit a 
physical behavior, i.e. they decrease. In the right plot of fig. 4, we display the differential 
cross sections for p* = 75 GeV and for several values of X. The symmetric-cut regime is 
here reached when X — )• 1, and in fact, as the values of X increases, the differential cross 
sections develop a negative tails for high Ayj^- 

The two plots in fig. 5, and similar ones that can be drawn for other values of X 
and p^, can help in the not-easy task of establishing what cuts to apply in an analysis if 
one wants to compare the experimental data to the NLO calculation using the alternative 
cuts of eq. (3.2). In the left plot, the cross sections plotted for several slices in Ayfj, 
become negative for values of p^ less than 75 GeV, fixing then a lower value below which 
the theoretical distributions cannot be trusted. Similarly, the right plot, where the cross 
sections are plotted as a function of X at fixed value of p* = 75 GeV, gives an indication 
of the upper value of X above which the NLO predictions become unreliable, for different 
slices of A^/fb- 

^While in principle X can vary from to 1, it can be shown that, for a partonic 2 — >■ n process, the 
lower limit is not 0, but l/(n — 1), so that, in our study (n = 3), 1/2 < X < 1. 
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Figure 5: Cross sections for different slices in Ayn, as a function of p* a,t X — 2/3 (left plot) and 
as a function of X a,t — 75 GeV (right plot). 

4. Probing higher-order corrections: a comparison among POWHEG, HEJ and 
NLO results 

In the previous section, we investigated various sets of cuts for avoiding regions of the phase 
space where the fixed NLO calculation becomes unreliable. We are now in a position to 
compare two theoretical approaches to dijet production that include higher order effects: 
POWHEG and HEJ. These two approaches are clearly very different in their description of QCD 
radiation. Nevertheless, for several kinematic distributions (see for example ref. [22]), the 
predictions from POWHEG and HEJ are very similar. We investigate here a few observables 
which could expose the differences in the two approaches and we compare them with the 
fixed NLO results. In order to avoid biasing our event sample, we impose a minimal set of 
asymmetric cuts 

p^T>35GeV, p^^>45GeY, \yj\ < A.7 , (4.1) 

i.e. all jets are required to have a minimum transverse momentum of 35 GeV, and the 
hardest-jet transverse momentum, ) is required to be greater than 45 GeV. In order to 
comply with the experimental acceptance, all jets are further required to have an absolute 
rapidity \yj\ less than 4.7. Jets are defined according to the anti-fcT jet algorithm, with 
radius R = 0.5. Only events with at least two jets fulfilling eq. (4.1) are kept. We stress 
that neither the POWHEG nor the HEJ descriptions suffer from using symmetric cuts, since 
they do a partial resummation of the large logarithmic terms. However, in order to have a 
sensible fixed NLO cross section to compare with, we impose asymmetric cuts. 

In the following, we compare the cross sections computed with a fixed NLO calculation 
and with HEJ, with the results obtained analyzing 14M events generated by the POWHEG 
BOX, at the level of the first-emission and after the shower performed by PYTHIA ^. The 
renormalization and factorization scales have been chosen equal to the transverse momen- 
tum of the hardest jet in each event, for the HEJ predictions. For the NLO computation 
(and for computing the POWHEG B function), scales are set to the transverse momentum 
of the underlying-Born configuration, as in the previous section. Scale-uncertainty bands, 

^We used PYTHIA 6.4.25 with the AMBTl tune. 
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Figure 6: The average number of jets as a function of A?/fb (left plot) and of Hj. (right plot), as 
predicted by a fixed NLO calculation, by POWHEG first emission, by POWHEG+PYTHIA and by HEJ. The 
dotted red lines around the HEJ prediction and the green ones around the NLO result are obtained 
by varying the renormalization and factorization scales by a factor of two around their central value. 



obtained by varying these scales by a factor of two in each direction, are shown for the 
NLO and HEJ results. The scales entering in the evaluation of parton distribution functions 
and of the strong coupling in the POWHEG Sudakov form factor are instead evaluated with 
a scale equal to the transverse momentum of the POWHEG hardest emission [29, 30]. 

The statistical errors (due to the numerical integration of the differential cross sections) 
on the theoretical predictions for the ratios studied in this section have been computed with 
the standard propagation of errors for the ratio of two uncorrelated quantities. We expect 
these errors to be an overestimation of the true ones, since, in our case, numerator and 
denominator are strongly correlated. To have a better estimation of the statistical error, 
we have applied a different method to the results produced by HEJ: we have divided the 
7 X 10^ generated events into 13 non-overlapping samples of equal size, and investigated 
the variation in the determination of the physical distributions of interest, when using all 
286 possibilities of choosing groups of 10 sub-samples out of the 13 ^. The estimate of the 
statistical fluctuation is determined by requiring that 68% of the results from these samples 
lies within the upper and lower boundaries. The numbers thus obtained agree well with 
the estimate obtained by using the standard variance form the same samples. 

In fig. 6 we plot the average number of jets as a function of the rapidity difference A?/fb 
between the most forward and most backward of the jets fulfilling eq. (4.1), on the left, and 
as a function of Hr^ = YIjPt the right. The wide-angle resummation implemented in 
HEJ produces more hard jets than POWHEG and the fixed NLO calculation, as Ayf^ increases. 
Both the NLO and the first-emission POWHEG results have at most 3 jets, so that the average 
number of jets cannot exceed 3, and are in good agreement. Additional jets are instead 

^The numbers 10 and 13 were chosen such that a reasonable convergence was obtained in each group, 
while maintaining a sufficient number of ways to pick 10 sub-samples for the statistics on the variation to 
be trustworthy. 
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Figure 7: Ratio of tlie inclusive 3-jet rate over the inclusive 2-jet one, as a function of Aj/ft (left plot) 
and Ht2 (right plot). The dotted lines display the results obtained by varying the renormalization 
and factorization scales by a factor of two around the central value, for the HEJ (red) and NLO 
(green) predictions. 

produced by the PYTHIA shower, so that the average number of jets is increased by roughly 
20% with respect to the NLO one, for Ayfb ~ 7. For the same separation in rapidity, the 
HEJ prediction is 45% larger than the NLO result, with a chance to distinguish among 
the three approaches. This same distribution has been investigated also in the analysis by 
ATLAS [22], but, as discussed earlier, the cuts applied in that analysis enhance the effects 
of collinear emissions, which are treated to leading-logarithmic order in POWHEG and partly 
in HEJ. In fact, the cuts in eq. (4.1) do not impose a large hierarchy of transverse scales 
and therefore avoid inducing dominance by collinear-soft emissions. The scale variation 
(dotted lines around the NLO and HEJ curves) are modest, of the order of a few percent. 
As a final comment, we note that the prediction from HEJ was found to be very stable 
against the effects of further showering [41] (using cuts very similar to those in eq. (4.1)). 

The dependence of the average number of jets on H^^ (right plot) displays a different 
behaviour: here the showered events have, on average, more jets than HEJ and the NLO 
results, as the sum of the transverse momentum of all the final-state jets increases. It 
is interesting here to comment on the NLO result obtained with the factorization and 
renormalization scales set to p^^/2, half of the transverse momentum of the underlying- 
Born configuration, i.e. the upper green dotted line in the plot. In fact, this quantity 
displays an unphysical behaviour: the average number of jets is greater than 3 for Hr^ > 
270 GeV. We will comment on this after the discussion of fig. 7, that suffers from the same 
problem. 

In fact, in fig. 7, we plot the ratio of the 3-jet inclusive cross section over the 2- 
jet one, as a function of Ay^, (left plot) and as a function of 11^2 (right plot), where 
Ht2 = Pt +Pt is the sum of the transverse momenta of the two hardest jets in the event. 
The same comments made about the left plot of fig. 6 apply here: the BFKL-inspired 
resummation implemented by HEJ produces more hard jets than the resummation of the 
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parton shower of PDWHEG+PYTHIA, for large rapidity separation between the most forward 
and most backward jet. For this distribution, an experimental analysis should then be able 
to distinguish between the HEJ and the POWHEG+PYTHIA predictions. 

In the right plot, the NLO ratio for the 3-jet inclusive cross section over the 2-jet one, 
plotted as a function of Ht2^ becomes unphysical (i.e. it gets greater than 1) when the 
factorization and renormalization scales are set to p}^^/2, for higher values of Ht2. This 
result is linked with the same unreliable behavior of the NLO distribution in the right plot 
of fig. 6, and deserves an explanation. For ease of notation we introduce the following 
shortcuts 

^mcl _ "-"2 ~mcl _ ""3 2) 

where the lower index 2 or 3 indicates the number of jets, and the differential cross sections 
are inclusive with respect to the corresponding number of jets. Together with inclusive 
cross sections, we define the exclusive ones, that will be designated with the upper label 
"excl". At fixed NLO, in dijet production, we have 

;:iincl ;:iexcl ;:iincl ;r.excl i ~excl {a q\ 

(^3 =^3 ' (^2 =^2 +<^3 ' (4.3j 

and we can relate the average number of jets with the inclusive 3-jet over 2-jet ratio 



Q ~ excl 1 Q ~excl pj-excl ;;:^incl 

z (T9 -I- ■303 =2+ =2 + — 

i- (J3 U2 O 2 



(jets) = .^^^^^ =2 + ^=2 + ^ (4.4) 

' excl I excl mcl incl ^ ' 



so that the two unphysical behaviors of the right plots in figs. 6 and 7 are strictly connected. 
The ratio 

Xincl Xincl 
U o (Jo 



- , ~ , -.^ (4.5) 

—.incl —.excl I —.mcl ^ ' 

"2 "^2 "T ''S 

can become greater than 1 for particular kinematic configurations only if the exclusive 2-jet 
cross section becomes negative at those phase space points. This happens, in our plots, 
when we choose as factorization and renormalization scale, for values of Ht2 ^ 

270 GeV. We have explicitly checked that the same behavior is observed if one sets the 
scale to be the hardest transverse momentum of the NLO partonic kinematics, a scale that 
is generally used for this kind of process. 

In the left plot of fig. 8, we plotted 1T2'"'', and, as expected, it becomes unphysical 
when the scale chosen is p}^^/2 and for values of Ht2 ^ 270 GeV. The explanation of this 
can be again traced back to the large logarithmic terms related to symmetric cuts and to 
the increase of the value of as, now evaluated at a smaller scale (the role played by the 
factorization scale would be more difficult to disentangle, since it involves the behavior 
of the parton distribution functions too). The 2-jet exclusive cross section always gets 
contribution from the Born and the virtual terms, irrespective of the value of Ht2, and 
from the part of the real-emission cross section that, at those kinematic points, is clustered 
into a 2-jet configuration 

ar' = + + = al (/i) {B + (fi) [V + R2]} , (4.6) 

where the notation is self-explanatory and wc put in evidence the appearance of the strong 
coupling constant q;s(a*) evaluated at the scale Since B and R2 are necessarily positive. 
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Figure 8: Exclusive 2-jet cross section as a function of i?T2, for three values of the renormalization 
and factorization scale fj, {p™,p"^/2, 2p™}, and for different set of cuts on the minimum-jet p^- 



coming from the square of the respective matrix element, it is the virtual term V that 
drives a^"^ to negative values. In other words, the R2 term becomes smaller and smaller 
if compared to the absolute value of V. At high values of Hj.2 , most of the events have 
three jets, with the two hard jets with Pt ~ ^t2/2, because, most hkely, the transverse 
momentum of the third jet is just high enough to pass the cuts in eq. (4.1). Since the 
minimum for the jets is 35 GeV, the imbalance between the two hardest jets is small, 
and the situation is equivalent to imposing symmetric cuts. On the other hand, if we 
increase to 100 GeV the minimum transverse momentum to define a jet, as in the right 
plot of fig. 8, the imbalance between the two hardest jet is no longer small compared to their 
average transverse momentum, and the event kinematics stays away from the symmetric- 
cut configuration. In addition, we have explicitly checked that, by keeping the p-r cuts of 
eq. (4.1) but restricting the jet-rapidity range to \yj\ < 2.5, the critical 7/^2 value is moved 
to higher values, i.e. ~400 GeV, implying that the unphysical behaviour depends on both 
the transverse momentum and the rapidity cuts. 

Turning now to the predictions from HEJ and POWHEG+PYTHIA in the right plot of 
fig. 7, we can see that they are in remarkable good agreement, both close to the POWHEG 
first-emission result, implying that the first POWHEG emission has the strongest impact on 
this distribution, while the subsequent shower has a milder effect. Before leaving this 
discussion, we would like to point out that, when using the cuts reported by the ATLAS 
Collaboration [23] 

pi^ > 60 GeV, p^' > 80 GeV, \yj\ < 2.8, (4.7) 

no significant difference has been found between the NLO result and the other three curves, 
so that the comparison the ATLAS Collaboration did is meaningful. 

As a last example of kinematic distribution that displays different behaviour if evalu- 
ated at NLO or using POWHEG or HEJ, we plot, in fig. 9, the average value of cos(7r — (/ifb), 
where (pit, is the azimuthal angle between the most forward and backward jets, as a func- 
tion of their rapidity separation Ayn,. For dijet events at tree-level, (/>fb = vr, since the 
two jets are back-to-back, and the average value of the cosine is 1. Deviation from 1 then 
indicates the presence of additional emissions, so that this kinematic distribution carries 
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Figure 9: The average value of cos(7r — (pf^^) as a function of Aj/ft, wtiere is tlie azimutlial angle 
separation between the most forward and most backward jet. The dotted red and green lines are 
obtained by varying the renormalization and factorization scales by a factor of 2 in both directions 
around the central value. 



information on the decorrelation between the two jets. This observable has been promoted 
for a long time as a good discriminator between descriptions with and without a systematic 
evolution in rapidity. It has therefore also been studied with full detector simulation in 
ref. [44]. However, the striking prediction from pure leading logarithmic BFKL evolution 
of an azimuthal decorrelation much larger than that which is obtained in a parton shower 
or fixed order formulation has been brought into question for some time by the inclusion of 
subleading corrections [45, 46]. This quantity is more inclusive than the average number of 
jets as it is sensitive also to emissions below the jet p-r cut. The higher radiation activity in 
POWHEG+PYTHIA and in HEJ, with respect to the fixed NLO and the POWHEG first-emission 
results, is clearly visible in the figure: the stronger jet activity produced by HEJ at higher 
rapidity separation (see the left plot of fig. 6) lowers the average value of the cosine below 
the POWHEG+PYTHIA result. As expected, the average value predicted by the POWHEG first- 
emission and the NLO calculation is closer to 1, since they contain at most one radiated 
parton. 

5. Conclusions 

Recent analyses by the ATLAS and the CMS Collaborations of inclusive and exclusive dijet 
production showed a high level of agreement between the two very different approaches to 
the description of perturbative higher-order corrections implemented in the POWHEG BOX 
and in HEJ. 

Inspired by these results, we have discussed the results obtained using a fixed NLO 
calculation, HEJ and POWHEG+PYTHIA, in the description of several kinematic distributions. 
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selected in order to display more clearly the differences among the three approaches: the 
average number of jets, the ratio of the inclusive three-jet rate over the two-jet one, and 
the azimuthal decorrelation between the most forward and the most backward jet. 

In order to perform a meaningful comparison between the NLO predictions and the 
resummed results obtained by POWHEG and HE J, we first investigated the effect of a few sets 
of cuts, and wc discussed which of them do not undermine the convergence of the fixed NLO 
results. In addition to the commonly used asymmetric cuts, we also considered a slightly 
different variant of them (the BFKL-inspired cuts of eq. (3.2)). It would be interesting to 
further explore the usage of these different cuts both theoretically and experimentally. 

Then we studied several kinematic distributions where the predictions of POWHEG+PYTHIA 
and HEJ differ significantly. While the limitations of the NLO calculation are clearly visible 
when we probe regions of the phase space where multi-jet emissions become important, the 
predictions of POWHEG+PYTHIA and HEJ are distinguishable when dealing with the average 
number of jets and the ratio of the inclusive three-jet production over the inclusive two- 
jet production, when studied as a function of the rapidity separation of the most forward 
and the most backward jet. Less marked differences are found when these quantities are 
plotted as a function of the sum of the transverse momenta of all the jets, or as a function 
of the transverse momenta of the two leading jets. Contrary to these findings, the study 
of the azimuthal decorrelation of the most forward and backward jet turned out to be less 
promising in distinguishing the two descriptions given by POWHEG and HEJ. 

We hope that an experimental measurement of dijet data collected at the LHC will 
follow, in order to investigate the quality of the theoretical understanding of these kinematic 
distributions. 
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